A fast eikonal equation solver using the Schrödinger wave equation

نویسندگان

  • Karthik S. Gurumoorthy
  • Anand Rangarajan
چکیده

Abstract. We use a Schrödinger wave equation formalism to solve the eikonal equation. We show that a solution to the eikonal equation is obtained in the limit as ~ → 0 of the solution to the corresponding linear Schrödinger equation. The Schrödinger equation corresponding to the eikonal turns out to be a generalized, screened Poisson equation. Despite being linear, it does not have a closed-form solution for arbitrary forcing functions. We use a standard perturbation analysis approach to derive a new algorithm which is guaranteed to converge provided the forcing function is bounded and positive. The perturbation technique requires a sequence of discrete convolutions which can be performed in O(N logN) using the Fast Fourier Transform (FFT) where N is the number of grid points. A major advantage of our approach over most other methods is that we do not require a spatial discretization of gradient operators and this contributes to the increased accuracy of our technique. The solution to the eikonal solution is recovered from the exponent of the wave function. Since the wave function is computed for a small but non-zero ~, the obtained solution is an approximation. We provide evidence for the usefulness of our technique by comparing the results of our approach with those obtained from Hamilton-Jacobi solvers such as the fast sweeping algorithm and the Dijkstra single source shortest path algorithm as well as with exact solutions. The latter are available for the Euclidean distance function problem—a special case of the eikonal when the forcing function is set to one.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

A Schrödinger Wave Equation Approach to the Eikonal Equation: Application to Image Analysis

As Planck’s constant ~ (treated as a free parameter) tends to zero, the solution to the eikonal equation |∇S(X)| = f(X) can be increasingly closely approximated by the solution to the corresponding Schrödinger equation. When the forcing function f(X) is set to one, we get the Euclidean distance function problem. We show that the corresponding Schrödinger equation has a closed form solution whic...

متن کامل

Paraxial geometrical optics for quasi-P waves: theories and numerical methods

The quasi-P wave in anisotropic solids is of practical importance in obtaining maximal imaging resolution in seismic exploration. The geometrical optics term in the asymptotic expansion for the wave characterizes the high frequency part of the quasi-P wave by using two functions: a phase (traveltime) function satisfying an eikonal equation and an amplitude function satisfying a transport equati...

متن کامل

Comparison of parallel solution techniques for the Eikonal equation

There have been many approaches to this end. An electrical wave travelling through heterogeneous heart tissue is usually simulated using a sophisticated system of coupled ODEs and PDEs. These equations are solved on a realistic geometric mesh with numerical schemes (an example using finite element techniques is shown in §1.1). In cases where a fast approximation to the wave dynamics is needed, ...

متن کامل

Optimal Control of Light Propagation Governed by Eikonal Equation within Inhomogeneous Media Using Computational Adjoint Approach

A mathematical model is presented in the present study to control‎ ‎the light propagation in an inhomogeneous media‎. ‎The method is ‎based on the identification of the optimal materials distribution‎ ‎in the media such that the trajectories of light rays follow the‎ ‎desired path‎. ‎The problem is formulated as a distributed parameter ‎identification problem and it is solved by a numerical met...

متن کامل

A uniformly second order fast sweeping method for eikonal equations

A uniformly second order method with a local solver based on the piecewise linear discontinuous Galerkin formulation is introduced to solve the eikonal equation with Dirichlet boundary conditions. The method utilizes an interesting phenomenon, referred as the superconvergence phenomenon, that the numerical solution of monotone upwind schemes for the eikonal equation is first order accurate on b...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:
  • CoRR

دوره abs/1403.1937  شماره 

صفحات  -

تاریخ انتشار 2009